Influence of dephasing process on the quantum Hall effect and the spin Hall effect 
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We study the influence of the phase relaxation process on Hall resistance and spin Hall current 
of a mesoscopic two-dimensional (2D) four-terminal Hall cross-bar with or without Rashba spin- 
orbit interaction (SOI) in a perpendicular uniform magnetic field. We find that the plateaus of the 
Hall resistance with even number of edge states can survive for very strong phase relaxation when 
the system size becomes much longer than the phase coherence length. On the other hand, the 
odd integer Hall resistance plateaus arising from the SOI are easily destroyed by the weak phase 
relaxation during the competition between the magnetic field and the SOI which delocalize the edge 
states. In addition, we have also studied the transverse spin Hall current and found that it exhibits 
resonant behavior whenever the Fermi level crosses the Landau band of the system. The phase 
relaxation process weakens the resonant spin Hall current and enhances the non-resonant spin Hall 
current . 
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PACS numbers: 73.23.-b, 72.25.Dc, 73.43.-f, 71.70.Di 



I. INTRODUCTION 

When the 2D electron system is subjected to a strong 
perpendicular magnetic field, the energy spectrum be- 
comes a series of impurity broadened Landau bands with 
extended state in the center— £ of each Landau band and 
the localized state at the band edges. This gives rise 
to the integer quantum Hall effect (IQHE);^ in which 
the Hall conductance is quantized and jumps from one 
quantized value to another when the Fermi energy sweeps 
through the impurity broadened bands or the extended- 
states. Experimentally^ the quantized unit of Hall con- 
ductance h/2e 2 can even be specified in parts per million, 
which becomes a resistance standard that is insensitive 
to the particular sample and details of its fabrication, 
because of the spatial separation in the transport states 
with the opposite velocity and consequently the incred- 
ibly long mean free path in the quantum Hall sample. 
Since the IQHE is discovered in 1980, it has been exten- 
sively studied in the past several decades ,£iL£iiliI£ and the 
many characteristics of IQHE has been well understood 
now. 

In addition to the quantum Hall effect, the spin Hall 
effect (in which the longitudinal electronic field induces 
the transverse spin current), especially the intrinsic spin 
Hall effect (SHE) has recently been intensively studied. 
Different from the extrinsic SHE which is due to the 
spin dependent scattering^ the intrinsic SHE is origi- 
nated from spin-orbit interaction (SOI). It was predicted 
first by Murakami et.alr^- and Sinova et.alM- in a Lut- 
tinger spin-orbit coupled 3D p-doped semiconductor and 
a Rashba spin-orbit coupled 2D electron gas, respectively. 
Subsequently, many related works are focused on intrin- 
sic SHE^>l£»a i8 . 19 i 20 i 21 i 22 Sheng et alM investigated 
the SHE in the mesoscopic 2D junction with Rashba SOI, 
and found that the SHE can survive in the mesoscopic 
systems with weak disorder. Xing et a££ found that the 



intrinsic SHE is dominated by the extend states, which 
is different from the IQHE. Besides, the out-of-plane and 
in-plane component of transverse spin Hall current was 
studied in a ballistic 2D finite electron system . 17 ' 18 On 
the experimental side, Kato et al^ and Wunderlich et 
aij^i have observed the transverse opposite spin accumu- 
lations near two edges of their devices when the longitu- 
dinal voltage bias is added. In addition, the reversal SHE 
was observed by measuring an induced transverse voltage 
in a diffusive metallic conductor when a longitudinal net 
spin current flows through iti^ 2 - 

Although the IQHE is insensitive to the particular 
sample and details of its fabrication, it can be trans- 
formed to the insulating regime with the weak magnetic 
field or strong impurity scattering. In fact, global phase 
diagram of transitions between the quantum Hall states 
and the insulator (or localized) state has been studied for 
the quantum Hall system in the tight-binding model 2 ^ 
and in 2DEG models A phase diagram for the meso- 
scopic SHE has also been proposed by Qiao et a/.. 19 These 
works showed both IQHE and SHE can survive in weak 
disorders and IQHE is more robust than SHE in resist- 
ing impurity scatterers. However, in a realistic sample, 
there exists both impurities or rigid scatters that main- 
tains the phase coherence and the dynamic scatterers like 
lattice vibration (photons) and electron-electron interac- 
tions that induce the phase-relaxation (PR). Hence it is 
interesting to ask to what extent the IQHE or SHE can 
survive in the presence of PR processes? 

In this paper, we study influence of the PR processes 
on the IQHE and SHE based on non-equilibrium Green's 
function (NEGF). We consider a 2D mesoscopic device 
which is sketched in the inset of Fig. 1(c): The central 
square ballistic region is connected to the four ideal semi- 
infinite lead-1,2,3 and 4 with the width W. The whole 
system, including the central region and the four leads, 
lies in the x, y— plane. A magnetic field B z is applied in 
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the positive z— direction. The PR processes in the cen- 
tral region are phenomcnologically simulated by intro- 
ducing the virtual leads 2 - , through which the electrons 
lose their phase memory. This method was first intro- 
duced by Biittiker in 1986. Moreover, there was another 
method to mimic the dephasing process by Datta and his 
co-worker.— This method provides a NEGF-based phe- 
nomenological model that is comparable to the virtual 
leads method with conceptual and numerical simplicity. 
In the following we will use the virtual leads method (also 
known as Biittiker probe model in some references) to 
mimic the dephasing processes. For the system with- 
out SOI, the longitudinal current or conductance J\ (we 
have set the bias V\ — V3 = 1) in the lead-1, transverse 
Hall voltage Vh = V2 — V4 and the step-like Hall resis- 
tance ph = Vh/J\ are calculated numerically with the 
aid of Green's function method. The results show that 
IQHE can survive at strong PR process. In particular, 
the quantized plateaus of the Hall conductance can be 
kept well even when the PR process is so strong as to 
completely relax the transport current. In the presence 
of the SOI, spin degeneracy is broken and odd integer 
Hall plateaus emerge. These odd number edge states are 
easily destroyed by weak PR processes. In addition, we 
also investigate how the SHE is affected by the PR pro- 
cess. It is found that the spin Hall currents J2/4,s show 
the resonant behaviors when the filling factor changes 
from odd to even where the Fermi energy is in line with 
a branch of the eigen levels of the spin polarized system. 
Furthermore, PR processes weaken the resonant spin cur- 
rent but enhance non-resonant spin current. 

The rest of this paper is organized as follows: In Sec. 
II, the system Hamiltonian and the theoretical formula 
for calculating the Hall resistance and other quantities 
are presented. In Sec. Ill, we show the numerical results 
and some discussions. Finally, a brief summary is given 
in Sec. IV. 



II. HAMILTONIAN AND FORMULA 

Our system (in the absence of phase relaxation and 
disorders) can be described by the Hamiltonian H = 
(—ihV + eA) 2 /2m* + a(a x V y — <r y V x ) with e the elec- 
tron charge, m* the effective mass, a the strength of the 
Rashba SOI, and a x , y the Pauli matrix. Generally speak- 
ing, there are always impurities or disorders in the real- 
istic system giving rise to rigid scatterings which do not 
contribute to phase-relaxation process. This type of im- 
purity scattering has been extensively investigated using 
the model of the on-site white-noise potential V\ (i de- 
notes the lattice site) distributed uniformly from -D/2 to 
D/2M^&±2- On the other hand, in the presence of phase- 
randomizing collisions, such as the dynamic scattering 
by the lattice vibrations (phonon), the electron-electron 
scattering, and so on, the transport electron loses the 
phase memory due to these PR processes. Although 
a proper treatment of non-coherent transport requires 



advanced delicate concepts, the basic issues can be ac- 
counted for by introducing the virtual leads 2 ^ attached to 
the site i to mimic the phase-breaking process occurring 
at the site ii 2 ^ Then the tunneling electrons can escape 
from the site i into the virtual leads where the electrons 
lose phase memory completely and finally return back 
to the site i. We assume that the PR processes occur 
only in the central region and the leads are treated as 
the measurement terminal which is ideal and clean. 

For the central and four real leads, we introduce the 
tight-binding representation, and the virtual leads are 
assumed in the free-electron form and expressed in the 
A;— space, which is not necessary but for simplicity. Then, 
the Hamiltonian is written in the following form : 16 ' 19 

i 

+ a\(tao + iVRU x )ai+s y + H.c 
+ y^Q; [B z g s fi B <Jz/2}ai 

i 

+ [ £ fe a ife a ife + (^ a i a ife + H - c ) W 

i,k 

where the first term describes the nearest neighbor cou- 
pling and the Rashba SOI in the central region and real 
leads, in which, i = (i x ,i y ) describes the site of the 2D 
region shown in the inset of the Fig. 1(c). ai — [<Zi,-|-, Q,i,x] T 
is the annihilation operator of electrons on the lattice site 
i, and Sx, Sy are unit vectors along the x and y directions. 
a x,y,z are Pauli matrices and tro is a 2 x 2 unit matrix. 
t = ti 2 /2m*a 2 is the nearest neighbor hopping matrix el- 
ement with the lattice constant a. tf> = qB z a 2 jh is the ex- 
tra phase unit originated from the vector potential A and 
m comes from \ y = ma. In the presence of the magnetic 
field B = (0,0, B z ), we introduce the vector potential 
A = (— B z y, 0,0) in the Landau Gauge, so that the ex- 
tra phase <p occurs in the form e - ™^ in the x— direction. 
Vr = ce/2a denotes the Rashba SOI strength in the tight- 
binding representation. In our model, the magnetic field 
B z is uniform in the whole system, including the four 
leads and the center region. While the SOI strength Vr 
exists only in the longitudinal lead-1, 3 and the central re- 
gion, and Vr is set to zero in the transverse lead-2,4. The 
second term in Eq.(l) denotes the Zeeman split where 
B Zl g s and fiR are the magnetic field along z— direction, 
the Lande g factor, and Bohr magneton, respectively. Fi- 
nally, the last term in Eq.(l) represents the Hamiltonian 
of the virtual leads (described in the k— space) and their 
coupling to the central site i. = [aifcj , aik,i] T is the 
annihilation operator of the electrons in the virtual leads, 
where i signs the positions of central region. Every cen- 
tral site is coupled by a virtual lead, so there totally are 
N = W x W virtual leads, with N being the site num- 
ber in the central region and W being the width of the 
central region. 

The charge current in real leads J r (r = 1, 2, 3, 4) and 
in virtual leads J v (v — 1,2, ...N) can be obtained from 
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the Landauer-Biittiker formula: 27 



j p = jY, t ™( v p- v ^> 

9#P 



(2) 



where p, q £ r or v, V p is the bias in the lead-p and 
T Pi9 is the transmission coefficient from the lead-g to the 
lead-p. The transmission coefficient can be calculated 
from T pq = Tr[r p G r r,G°], where the line-width func- 
tion F p = — £pt) with Sp the retarded self-energy 
and the Green's function G r = [G a ]t = {E F I - H - 
^2 p where I is the unit matrix with the same di- 

mension as that of H . In addition, in order to inves- 
tigate the SHE, we also need to calculate the spin Hall 
current J p jS in the transverse lead-2 and 4 for the system 
with SOI. Because we have set Vr = in lead-2 and 4 so 
that a is a good quantum number, the particle current 
J pa in the lead-p (p — 2, 4) with spin index a (er ="\ or 
j.) can be obtained from the Landauer-Biittiker formula: 



v q ). 



(3) 
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FIG. 1: (Color online) The total current Jt, the coherent 
component J c , and the non-coherent component J nc vs. the 
PR strength Ta in the two-terminal non-SOI system for the 
magnetic field B — 0.5 (panel (a)) and B = 0.1 (panel (b)) 
and the system size W = 20 (green or gray lines) and W = 40a 
(black lines). In the Panel (c), the PR strength Yd vs. phase 
coherence length L^, is plotted. Inset of panel (c): Schematic 
diagram for the mesoscopic four-terminal device. 



The quantities in Eq.Q are the same as that in Eq.([2]), 
except that here T paq is the transmission coefficient 
from the lead-g to the lead-p with spin er, and T pa ,q = 
Tr[T ptT G r T q G a ] withr p(T = i(S;„-S;t). After obtain- 
ing the particle current J pcr , the spin current J PjS can be 
easily obtained as J p , s = (h/2)[J p ^ — J p jJ. 

In our calculation, the external bias is applied in the 
longitudinal leads-1 and 3 with V\ = 0.5 and V3 = —0.5, 
thus electrons obtain a velocity v x along the x— direction. 
With the perpendicular magnetic field, the tunneling 
electrons deflect to the transverse direction (y-direction) 
and generates the charge pile up at the vicinity of the 
lead-2 and lead-4. It consequently leads to opposite 
transverse Hall voltage V2 — — V4 in the lead-2 and lead- 
4 which is calculated by requiring the boundary condi- 
tion J2 = J4 = since the lead-2 and 4 only act as 
the voltage probes. Furthermore, the electrons can only 
lose their phase memory by escaping into or coming back 
from the virtual leads, and do not contribute net cur- 
rent to the virtual leads, so there are N extra boundary 
conditions J v = with v = 1,2, ...N. With the trans- 
mission coefficient T pq , the input parameters V\ and V3, 
and the boundary condition J P= 2A, V — 0, we can get the 
transverse bias V24 (V2 = — V4) in the lead-2, 4 and the 
longitudinal current J1.3 (Ji = — J3) in the lead- 1,3 us- 
ing Eq.([2]). Consequently the Hall voltage Vh = V2 — V4 
and Hall resistance pn = Vh / J\ arc obtained straightfor- 
wardly. On the other hand, when investigating the SHE 
we use the boundary condition V% = V4 = instead of 
J2 = J4 = 0. Then, the particle current 3 p . a and conse- 
quently the spin current J p . s in the lead-2, 4 can also be 
calculated easily. 



III. NUMERICAL RESULTS AND DISCUSSION 

In the numerical calculation, we fix the Fermi energy^ 
Ep = —3t which is near the band bottom — 4t, and take 
t as energy unit. Here t — h 2 /(2m*a 2 ) is about 5meV 
while taking m* = 0.05m e (the mass of free electron) and 
the lattice constant a — 12.5nm. The size of finite sample 
W is confined to W = 40a and 20a. In the presence of 
magnetic field B Zl there is an extra phase unit <p = j^Ba 2 
with the vector potential A = (— B z y, 0, 0). When <f> — 
1, B = h/(ea 2 ), so we take h/(ea 2 ) as the unit of the 
magnetic field B, which corresponds to B — 4.2T and 
HbB = 0.05i. We set Lande g factor g s = 2. Thus 
Zeeman splitting }^9s^bB — 0.05i. Moreover, in order to 
generate the edge states, the cyclotron radius r c = v/u) c 
must satisfy r c < W/2, then the magnetic field 2/W < 
B is needed. In the SOI system, Vr = xt corresponds 
to the strength of Rashba SOI a — 2a Vr ~ 1.25a; x 
lQ~ w eVm and spin precession length (over which the 
precessing angle n is generated) Lso — 7ra/(2x). Finally, 
as a check for our computer program, we have calculated 
the case in which the magnetic field and the PR process is 
absent but the spin-orbit interaction is present, the same 
result as in the Fig.l of the first reference in Ref.14 can 
been obtained. 

In the experiment, the phase coherent length is a 
observable parameter and is used to describe all kinds of 
dephasing processes. Under different experimental con- 
ditions, there arc different dephasing processes and cor- 
respondingly different L$. So in the following, we will 
study the relation between the PR strength and the 
phase coherent length L$. Here we define as the 
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length, through which the transport electron has 50% 
probability to lose its phase memory and 50% probabil- 
ity to keep the phase coherence. In the presence of PR 
processes T^, the current generally consists of both phase 
coherent part and phase incoherent part. In order to es- 
timate the phase coherence length L§ for a given PR 
strength Yd, we consider a two-terminal structure (de- 
coupled with the lead-2,4), in which the electrons can di- 
rectly flow from lead-1 to lead-3 or indirectly from lead-1 
to lead-3 through virtual leads, the former contributes to 
the phase coherent current and the latter the incoherent 
part. Fig. la and b show the coherent component, the 
incoherent component, and the total current versus 
for the different system sizes (i.e. the system lengths) 
W = 20a and W — 40a. The total current decreases 
slightly with the increasing Yd, since the edge states car- 
rying the transport electrons are slightly destroyed to 
the extent proportional to Yd- At = 0, the incoherent 
component is zero and the coherent component is equal 
to the total current. With the increase of Yd, the inco- 
herent component increases and the coherent component 
decreases. At a certain Td, i.e., at the crossing point of 
the two curves, the incoherent component is just equal to 
the coherent component. This means that for this crit- 
ical Tdc the system length W is just equal to the phase 
relaxation length L^. Therefore the relation between L$ 
and the PR strength Yd can be obtained. Fig.lc shows 
Yd versus L$ for different magnetic fields B. The phase 
coherence length increases monotonically when Yd de- 
creases. Obviously, when Yd — > 0, — > oo. When Yd is 
not very small (e.g. Yd > 0.05t), the bigger the magnetic 
field B, the longer is. This is because for the stronger 
B the edge states are more robust against the PR process. 
In particular, for the given system size W = 40a and 20a 
used in our following calculation, the critical Yd c ~ 0.02 
and 0.05, respectively, beyond which the phase coherence 
length is smaller than the system size. 



A. the IQHE 

In this subsection, we study the influence of the PR 
processes on the IQHE. Firstly, we consider the system 
without Rashba SOI. In the absence of PR processes 
(r<j = 0), all of the electrons are carried by the edge 
states, the electrons traverse clockwisely along the sam- 
ple edge due to the Lorentz force qv x B (here B = B z 
is along the positive z-direction), and the only nonzero 
matrix elements of the transmission coefficient matrix T 
are T2i,T32,T43 and T14 with integer values. When PR 
processes exist (Yd ^ 0), the edge state is partially de- 
stroyed. As a result, the elements ?2i,T32,T43 and T14 
of the transmission coefficients, which denote the edge 
states, deviate from integer values, and the other ele- 
ments such as T\2 become nonzero. In Fig. 2, the current 
J 1 in the lead-1, the transverse Hall voltage Vh and Hall 
resistivity pn or 1/ph versus magnetic field B or X/B are 
plotted for the different r^. From Fig. 2, we can see when 
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FIG. 2: (Color online) Panel (a)-(c): the current in the lead-1 
Ji , the transverse Hall voltage Vh and Hall resistivity pn vs. 
magnetic field B for the different PR strength Yd, with Vr = 
and W = 40a. Panel (d) magnifies the marked region in the 
panel (c). Inverse of Hall resistivity X/ph vs. the inverse of 
the magnetic field X/B in the non-SOI system (Vr = 0) for 
the different system size W = 40a (panel (e)) and W = 20a 
(panel (f)). 



the magnetic field (B < 0.05) is too small to form edge 
states, Ji changes slowly for the variable B and decreases 
rapidly with Yd, the Hall voltage Vh and Hall resistance 
Ph increase linearly with B, which is in agreement with 
the results of semiclassical Drude model. In the following, 
we focus on the high fields B case, i.e., B is large enough 
to separate the flows with the opposite velocity^ and 
the system is in the quantum Hall regime. When = 0, 
the Hall voltage Vh = V% — V4 = 1, J\ and consequently 
Ph exhibit plateaus. This can be understood using the 
well known picture of Landau level (LL) and the edge 
state 30 . In the presence of PR process (r^ 7^ 0), J\ se- 
riously deviates from the even integer plateaus and Vh 
is no longer a constant. However pn hardly changes and 
still keeps plateaus even when Yd is much bigger than the 
critical value Ydc (see Fig. 2c and d). This means that al- 
though the system is in strong PR regime, the IQHE can 
survive and is rather robust against the PR processes. 
From Fig. 2(d), we can see that the the plateaus are first 
destroyed at the band edges for large Yd- This is because 
at the band edges of plateaus, Ep is closer to LLs than 
at the band centers of plateaus. As a result, the electrons 
are easier to be relaxed to LLs leading to a smaller pn ■ In 
Fig.2e and f, we plot X/pn versus X/B for very large Yd- 
The results show that it is more difficult to destroy the 
plateau at larger magnetic field B or larger sample size. 
So those plateaus can survive at a bigger Yd, i.e., they 
have stronger ability to resist the PR processes. This is 
because the energy spacing AE between the nearest LLs 
is larger for the larger magnetic field B, and the edge 
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FIG. 3: (Color online) Panel (a): inverse of Hall resistivity 
1/pH vs. the inverse of magnetic field 1/B for the different 
Rashba SOI strength Vr at F d = 0. Panel (b): l/p H vs. 1/B 
with Vr — 0.05t and 0.5t for different P^. The system size 
W = 40a. 



FIG. 4: (Color online) Panel (a): The transverse spin current 
J2,s vs. inverse of magnetic field 1/B for the different Rashba 
SOI strength Vr. Panel (b): The transverse spin current 
J2, s vs. the Rashba SOI strength Vr near the resonant peak 
1/B ~ 5.f in the panel (a). The other parameters are r d = 
and W = 40a. 



states carrying the opposite current are separated at a 
larger distance for the larger sample. In particular, the 
first plateau at 1/B < 3 can keep well even when Td 
reaches 10 which is two orders of magnitude larger than 
the critical value Tdc = 0.02 or 0.05. 

Next, we study the system with Rashba SOI Vr. In 
Fig. 3a, 1/pH versus 1/B for different Vr in the absence 
of PR processes (r^ = 0) is plotted. It is interesting that 
although the SOI induces the extended states^ 1/ ' pu 
is still quantized and the integer quantum plateaus still 
remain. At Vr — 0, there are only even quanta for 1/ Ph- 
However, as long as Vr ^ 0, odd quanta emerge because 
the spin degeneracy is destroyed. The width of the odd 
plateau is wider for the larger Rashba SOI strength Vr 
or for the smaller magnetic field B (see Fig. 3a). While 
at the large Vr, the width of the odd plateau can be in 
same order with, or even wider than, the width of the 
even plateau. 

In the following, let us study the influence of the PR 
processes on the plateaus of the IQHE with non-zero 
Rahsba SOI (Vr ^ 0). Fig.3 shows l/p H versus 1/B for 
different SOI strengths Vr and PR strengths Td- When 
Vr is small, for example Vr = 0.05i, the even plateaus 
originated from the edge states are still present. The 
odd plateaus due to SOI, however, are quickly washed 
out. While for the strong SOI (e.g. Vr = 0.5t) case, the 
SOI that favors the extend states dominates the mag- 
netic field that favors the edge states, both even and odd 
plateaus are destroyed by the weak PR process. Notice 
that for the non-SOI (Vr = 0) case the even plateaus can 
be kept at very large Td which can even reach 10 as shown 
in Fig.2e and f. The fact that the even plateaus in the 
SOI system are not as robust against the PR processes as 



that in the non-SOI system, indicates that the SOI weak- 
ens the ability of resisting the PR processes. In fact, the 
Rashba SOI coefficient a is usually less than 10 _11 eVm 
for a general 2D electron gas, and the corresponding Vr 
is less than 0.08i. With this Vr (e.g. Vr = 0.05i), the 
even plateaus can still survive at Td = 0.1 which is much 
larger than the critical value Td c - 



B. the SHE 

In the system with SOI, the SHE occurs, in which 
a pure and non-dissipating transverse spin current can 
be generated when a longitudinal electric field or bias 
is applied. Recently, the SHE has been extensively in- 
vestigated by a great deal of works as mentioned in the 
introduction-ii^^iii^^^iSi^ Here we mainly study 
how the SHE is affected by the magnetic field B, and, 
in particularly, the PR processes r^. Firstly, the case of 
= is studied. The transverse spin current J2.S — 
— J^s versus 1/B for different Vr is plotted in Fig. 4(a). 
An interesting feature is that the spin current J2.S shows 
a resonant behavior, when the quantized 1/ph changes 
from the odd plateau to the even plateau where the Fermi 
level is in line with the one of energy eigenvalues of the 
spin degenerated system. The origin of the resonant spin 
current will be discussed at last paragraph in this sec- 
tion (see Fig. 6). The spin Hall current J 2jS is quite large 
when 1/B is near the resonance, but is very small when 
1/B is far away from it. For Vr = 0.2, the resonant po- 
sitions are about 1/B — 3.1t, 5.1i, • • • , and for Vr = 0.1, 
they are about at 1/B = 3.0i, 5.0t, 9M, 7.0t, ■ ■ ■ . In 
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FIG. 5: (Color online) Panel (a): the transverse spin current 
J2, s vs. inverse of magnetic field 1/B for the different PR 
strength Td at Vr = 0.1. Panel (b) and (c): the transverse 
spin current J2,s vs. the Rashba SOI strength Vr for the 
different PR strength Td at the off-resonant peak 1/B = 2 
(b) and resonant peak 1/B = 3.1 (c) in the panel (a). The 
system size W = 40a. 



Fig. 4(b), we fix the magnetic field 1/B at a resonant 
point (1/B — 5.1) and plot J2.S versus Vr. The results 
show that J2, s is randomly distributed at the large Vr 
because the extended states are dominant. On the other 
hand, when Vr < 0.5i, J2, s is regular. For 1/B = 5.1, 
J2 : s is resonant for Vr = 0.2t, and there is a stable in- 
terval near Vr = 0.2t (see the red dash dotted line in 
Fig. 4b). When 1/B deviates from the resonant point 
5.1, J2, s is rapidly decay from the resonant Vr (see black 
dotted line and green solid line in Fig. 4b). 

In the following, we study the influence of the PR pro- 
cesses Td on the transverse spin current J2. s - In Fig. 5(a), 
we plot J 2, a vs the magnetic field 1/B for different Td- 
It shows that the PR process Td suppresses the spin cur- 
rent J2, s when 1/B is near the resonant points, but en- 
hances J2,s in the off-resonant region. Next, picking up 
the off-resonant position (1/B = 2) and resonant position 
(1/B = 3.1), versus the SOI strength Vr for different 
Td are plotted in Fig. 5b and c, respectively. Similar to 
Fig. 5 (a), the off-resonant spin current is enhanced by the 
PR processes (see Fig. 5b). On the other hand, when the 
magnetic field is at the resonant point (e.g. 1/B = 3.1 
in Fig. 5c), the resonant spin current with Vr from 0.2 
to 0.4 is suppressed (see Fig. 5c). Moreover, when Vr is 
very large (e.g. Vr > 0.5), the spin current enters the 
chaotic regime and it depends on the PR processes Td in 
a random fashion. 

Finally, we also study the influence of the Zeeman ef- 
fect on the SHE. Fig. 6(a) and (b) show the spin current 
J2,s versus the inverse of magnetic field 1 /B with the SOI 
strength Vr — and O.lt, respectively. In each panel, 
we also plot the corresponding eigenvalues (LLs) of the 
central scattering region (without leads) versus 1/B and 
assign the resonant positions (see the red dotted lines). 
When Vr — 0, there is only one group of peaks in the 
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FIG. 6: (color on line) The transverse spin current J2,s, the 
system Landau levels (LLs) vs. inverse of magnetic field 1/B 
while in the presence of the Zeeman effect, with the Rashba 
SOI strength Vr = (a) and Vr = O.lt (b). The other 
parameters are Td — and W = 40a. The Fermi level is fixed 
at E = -3t. 



curve of J2, s vs 1/B due to the Zeeman effect. While for 
Vr 7^ 0, the peaks appear in two groups corresponding 
to the Zeeman peak and resonant peak of spin current. 
These can be understood as follows. When considering 
the Zeeman effect, the density of state of electrons for the 
spin-up and spin-down is different, so LLs are split into 
the spin-up and spin-down channels. For a given Fermi 
energy Ep, the number of LLs below Ep for the spin-up 
and spin-down states can be different. If the number of 
LLs below Ep are even, the spin current is zero. How- 
ever, if the number of LLs below Ep is odd, i.e., when Ep 
is between the split of the LLs, then the spin current is 
one half and the (Zeeman) peak emerges in the spin cur- 
rent. The positions of peaks are just the positions of LLs 
(see the red dotted line in fig. 6(a)) since the Zeeman split 
is very small (0.05t). On the other hand, in the presence 
of SOI with Vr 7^ 0, the LLs for up and down spin chan- 
nels in the presence of Zeeman term are mixed together 
and become two new LLs with different spin polariza- 
tions. One of LLs is strongly spin polarized while there 
is no spin polarization for the other LLs. These results 
are obtained by solving the Schrodinger equation with 
Vr ^ 0. As a result, the resonant spin current emerges 
when the Fermi level is just in line with the spin-polarized 
LL (see the red dotted line in Fig. 6(b)). So, except for the 
Zeeman peaks, there are another group of peaks which 
is originated from different physics in the non-zero SOI 
system. From Fig. 6, we can also see that the intervals 
of Zeeman peaks are unchanged, while the intervals of 
the resonant peaks are closer with the decreasing Vr (no 
shown), which gives an extra evidence that they come 
from the different physics. 
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IV. CONCLUSIONS 

The effect of PR process Td is investigated in the 2D 
four-terminal system with or without SOI. Without the 
SOI, the plateaus of the Hall resistance pn are found 
to be well kept even when the PR strength is very 
strong (i.e. the phase coherence length is much shorter 
than the size of system) . This means that the IQHE has 
very strong ability to resist the PR process IV Further- 
more, for the larger sample or the stronger magnetic field, 
the resistive ability of the PR processes Td is stronger. 
On the other hand, for the system with SOI, the odd 
integer plateaus of 1/pn are also appear. The odd in- 
teger plateaus due to the SOI can be destroyed even for 
the very weak r^, but the even integer plateaus can still 
survive in quite strong PR process Yd- Next, the SHE, 
i.e. the transverse spin current, is also studied in the 



system with SOI. It is found that the transverse spin 
current reaches the resonant pole when the Fermi level 
is just consistent with the one of the two branches of the 
energy eigenvalues of the system with SOI. The PR pro- 
cess weakens the resonant spin Hall current and enhances 
the non-resonant spin Hall current. In addition, we also 
study the properties of system with the Zeeman effect 
and find there are two group of peaks of the spin current 
originated from different physics. 
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